High quality asynchronous heralded single photon source at telecom wavelength 
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We report on the experimental realization and characterization of an asynchronous heralded 
single photon source based on spontaneous parametric down conversion. Photons at 1550 nm are 
heralded as being inside a single-mode fiber with more than 60% probability, and the multi-photon 
emission probability is reduced by a factor up to more than 500 compared to Poissonian light sources. 
These figures of merit, together with the choice of telecom wavelength for the heralded photons are 
compatible with practical applications needing very efficient and robust single photon sources. 
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With the present development of quantum communi- 
cation and computing technologies, including quantum 
key distribution pj and quantum teleportation the 
interest for true single photon sources is rising. Many 
different implementations have been investigated, using 
single molecule or atom excitation B, 1 3, B- B, color cen- 
ters in diamonds 0,0, quantum dots |9l I lOifllL or 
pulsed parametric down conversion sources |l4l |l5j . All 
theses solutions have various advantages and tradeoffs 
between high purity and efficient single photon produc- 
tion, repetition rate, wavelength of the photons, and ease 
of use. The aim of this paper is to show that a sponta- 
neous parametric down conversion (SPDC) source made 
of a bulk non-linear cristal at room temperature and a 
simple basic optical setup can be used to herald single 
photons at telecom wavelength in a very efficient way 
(see figure The photons are also directly available 
in a standard single-mode telecom optical fiber, making 
this source a good choice for quantum communication 
applications such as scalable quantum networks. In this 
context, the term "heralded" means that photons are not 
generated on demand, but instead an electric signal an- 
nounces the presence of a photon in a fiber. Indeed, as 
photons are created in pairs, the detection of one photon 
can be used to announce the presence of the complemen- 
tary photon |l6| . 

The source presented in this paper is an asynchronous 
heralded single photon source (AHSPS) because the 
heralding signals are not synchronized with a periodic 
clock, the SPDC pump being continuous. It exhibits a 
very good probability of producing one photon and a very 
low probability of producing more than one photon per 
heralding signal. The latter probability depends on the 
pump power applied to the crystal, and several measure- 
ments are presented to characterize this dependency. 
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FIG. 1: Schematic of the asynchronous heralded single photon 
source (DM: dichroic mirror, DF: neutral density filter for 
pump attenuation, LPF: low-pass pump filters, VA: variable 
optical fiber attenuator) 



Our AHSPS is made of two main parts, as depicted in 
figure The first consists in the SPDC photon pair cre- 
ation stage which consists in a type I KNb03 bulk non- 
linear crystal, pumped with a 532 nm continuous-wave 
laser (with maximum power of about 60mW). On one 
hand we require our source to produce single photons 
for long range telecom applications, and on the other 
hand we need the complementary photons to be pas- 
sively and efficiently detected for heralding purpose. The 
phase matching condition is therefore chosen to produce 
1550-810 nm photons pairs, with the photons at 1550 mn 
featuring a 6.9 nm spectral width. A dichroic mirror is 
used to separate the 810 nm photons from the 1550 nm 
photons into two distinct beams, enabling one to couple 
these two wavelength using separated optimized optics 
into single-mode fibers. More details about the core of 
this source can be found in The 1550 nm single 

mode fiber is the heralded single photons output port. 

The second part of the AHSPS is made of a actively 
quenched silicon avalanche photodiode (APD) (EG&G) 
connected to the 810 nm single-mode fiber. The TTL de- 
tection signal from the silicon counter is the asynchronous 
heralding electric signal. For characterization, the opti- 
cal variable attenuator in front of the silicon APD can be 
used to reduce the heralding signal rate without changing 
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the heralded photons statistics. For the experiment, the 
pump power is varied by inserting neutral density filters 
in the path of the pump laser beam. 

The heralded photon output port is connected to a test 
bench consisting of a gated detection Hanbury-Brown & 
Twiss type setup |l9j . In order to characterize the source, 
the relevant parameters P(l), P(2) and g^ 2 '(0) have in- 
deed to be measured. P(l) is the probability of hav- 
ing exactly 1 photon (and P(2) the probability of hav- 
ing more than 1 photons) in the photon output port per 
heralding signal and during a time corresponding to the 
detection gate width. Note that the Hanbury-Brown & 
Twiss setup does not enable to discriminate between 2 
and more photons events. By convention, we thus have 
P(0) + P(l) + P(2) = 1, where P(0) is the probability 
of having no photon at all per heralding signal. Finally, 
(0) is the standard second-order autocorrelation func- 
tion at time |2(| . 

The setup is made of a 50/50 fiber optical beam splitter 
and two InGaAs APDs operated in gated mode (idQuan- 
tique). Both of these detectors are triggered with the 
heralding signal of the source through a TTL fan-out. 
The effective gate width is about 2.5 ns and is chosen 
with respect to the jitter of Si APD, which is about 1 ns. 
A narrower gate width would not allow to detect all the 
heralded photons and would thus bias down the measured 
value for P(l) and P(2). The measured dark count prob- 
ability per gate is p A c = 35.1 ± 0.7 x 10~ 6 for detector A 
and p B ° = 7.4±0.3 x 10~ 6 for detector B. The dark counts 
of the AHSPS silicon detector (about 100 Hz) artificially 
increase the probability of having no photons in the out- 
put port, but this effect is measured to be negligible in 
the present case. The detection states of both detectors 
are acquired using a computer by means of a digital in- 
put pc card (PCI-6533 DIO-32HS, National Instrument) 
and a labview software (see figure EJ The triggering rate 
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FIG. 2: Schematic diagram of the test bench. The beam split- 
ter (BS) is considered as having a perfect 50/50 separation 
ratio and the fibers from it to the detectors are considered as 
lossless. This is possible because detectors efficiency are cali- 
brated including losses in the beam splitter and in the fibers 
(see text). 



of the APDs is kept under 50 kHz using the optical vari- 
able attenuator of the source. Note that this does not 
affect the quality of the source, but has to do with the 
test bench we use. Indeed, this low triggering rate and 
a 10/is external deadtime are used on both the InGaAs 
APDs to work in a regime where after-pulses counts are 



negligible. Because of these deadtimes, the heralding sig- 
nal of the source cannot be directly used to trigger the 
APDs state acquisition, but one must discard the cases 
where one or both detectors are disabled when a new 
trigger arrives. If not, an artificially lowered coincidence 
detection rate would be measured, and the quality of the 
source would be overestimated. This is done by using 
a signal named "gate out" on figure [3 which is synchro- 
nized with the effective application of bias voltage gate 
on a given APD, meaning that it is ready to detect her- 
alded photons. These signals are combined with a logic 
AND gate to form the acquisition triggering signal. Upon 
reception of a trigger signal, the detection states (0 or f ) 
of the two APDs are acquired independently. The raw 
output is a binary file containing the detection state for 
each heralding signal effectively taken into account. 

Using attenuated laser pulses, the efficiencies of both 
detectors under these particular working conditions are 
calibrated including the losses in the beam splitter and 
in the fibers connecting the APDs. It is thus possible to 
consider that the beam splitter is a perfect 50/50 separa- 
tor and that the fibers are lossless. The overall detection 
efficiencies are measured to be Pa = 8.4 ± 0.1% for de- 
tector A and (3 B = 9.6 ± 0.1% for detector B. 

For the present setup, one can show that the g^(0) 
parameter is equal to 



5 (2) (0) 



Pab 



_ 2P(2) 
Pa Pb ~ P(f) 2 



(1) 



where pab is the experimental probability of having si- 
multaneous detections in both detectors and pa, Pb the 
probabilities of having a single detection in the given de- 
tector A or B. Note that, as expected, this quantity does 
not depend explicitly on the detection efficiency. Refer 
to the appendix for details on the rightmost part of the 
equality and for details on obtaining P(l) and P(2) 
from the experimental data, pa and pb are obtained 
from the raw binary files by counting the number of reg- 
istered detection for the given detector {Na, Nb) divided 
by the total number of acquisition triggers (N t ). pab 
is obtained by counting how many times both detectors 
have registered a simultaneous detection (Nab) divided 
by the total number of acquisition triggers. 

For the measurement the SPDC source is first aligned 
without neutral density pump filters, by optimizing the 
collection of the 810 nm photons. A second alignment 
is done by maximizing the ratio of detection at 1550 nm 
over detection at 810 nm through the adjustment of both 
1550 nm and 810 nm photons collection. The highest ra- 
tio is usually found when the 810 nm detection rate is 
about 90% of the maximal rate. Different neutral den- 
sity pump filters are introduced one after the other in the 
path of the pump beam. For each filter, the pump power 
is measured just before the crystal, and the heralding 
signal rate Rh is recorded. When necessary, the vari- 
able attenuator is set to a non-zero value to keep this 
rate under 50 kHz for the acquisition. The results are 
summarized in table [I] In this table, P p is the pump 
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TABLE I: Experimental AHSPS parameters obtained with a detection gate width of 2.5 ns 



Pp 


Rh 


N t 


T 


N A 


N AB 


P(l) 


P(2) 


5 (2) (0) 




49 mW 


845 kHz 


6.0 x 10 7 


20 min. 


16 x 10 5 


1062 


0.59 


41 ± 1 x 10" 4 


22.6 ±0.5 x 10 


-3 


25 mW 


403 kHz 


2.8 x 10 7 


20 min. 


7.7 x 10 5 


255 


0.61 


19 ± 1 x 10~ 4 


10.2 ±0.5 x 10 


-3 


13 mW 


221 kHz 


1.6 x 10 7 


20 min. 


4.1 x 10 5 


82 


0.61 


10 ± 1 x 10~ 4 


5.3 ±0.5 x 10" 


3 


3.3 mW 


120 kHz 


1.7 x 10 7 


12 min. 


4.6 x 10 5 


78 


0.60 


8.0 ±0.9 x 10" 4 


4.4 ±0.5 x 10" 


3 


1.6 mW 


39 kHz 


1.6 x 10 7 


40 min. 


4.2 x 10 5 


35 


0.61 


2.5 ±0.6 x 10~ 4 


1.4 ±0.3 x 10" 


3 


0.9 mW 


18 kHz 


1.8 x 10 7 


85 min. 


4.4 x 10 5 


46 


0.58 


3.7 ±0.6 x 10 -4 


2.2 ± 0.4 x 10" 
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power and T is the experimental duration of the acqui- 
sition. -P(l), P(2) and g^(0) are net (dark count sub- 
stracted values) calculated from formulas developed in 
the appendix. These values are obtained from the av- 
erages of the slightly different results calculated using 
Na and Nb- As simultaneous detections are quite rare 
events, Nab is small. The statistical error on this quan- 
tity is the main source of error for P(2) and g( 2 \0). It 
can be relatively high but it is kept under an acceptable 
level by acquiring for relatively long times. The error on 
P(l) is negligible, because Na (and Nb) is large. The 
slight changes for P(l) values (less than 2% around the 
average value of 0.6) are explained by the fact that the 
neutral density pump filters slightly deviate the pump 
beam. This causes a change in the phase matching con- 
dition and consequently a decrease of the coupling ratio 
from its nominal value. As P(l) ^ -P(2), P(l) is to a 
very good approximation the heralded photons coupling 
ratio, which is sometimes called the heralding efficiency 
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FIG. 3: P(l) (left scale) and P(2) (right scale) displayed as 
function of the heralding rate. The solid lines are: average 
for P(l) and linear fit for P(2). As expected, P(2) is linearly 
dependent on heralding rate (and hence on the pump power) , 
while P(l) remains constant. 

These results should be compared to a Poissonian light 
source featuring the same value for P(l). The main ob- 
servation is that the multi-photon emission probability 



P(2) of our source is reduced by a factor l/g( 2 '(0) rang- 
ing from about 50 to more than 500, depending on the 
pump power. To the best of our knowledge, these are the 
highest factors reported so far. In fact, g( 2 '(0) and P(2) 
can be arbitrarily low for low pump powers, i.e. heralding 
rates, as P(2) is linearly dependent on the heralding rate 
(see figure|2Jl. Moreover, P(2) can be reduced by lowering 
P(l), i.e. misaligning the heralded photons coupling. A 
better figure of merit is therefore the ratio of the heralded 
photons rate over g( 2 \Q): 

P{l)R H = P{1?Rh 
.9(2) (0) 2P(2) 

Rh and P(2) are related to the physical parameters of 
the source as follows: Rh is the number of photon pairs 
created by second multiplied by the probability to detect 
the heralding photon. P(2) is the probability to find two 
photon pairs in a time interval At knowing that a trigger 
photon was detected, which is half of the probability that 
a photon pair is created during At, multiplied by the 
probability that two photons are coupled into the output 
fibre, P(l) 2 . We find 

R H = qHVHVsP P , P(2) = 0.5 P(l) 2 At Vs P p 

where qn is the coupling efficiency of the heralding pho- 
tons, t]h is the quantum efficiency of the heralding detec- 
tor, -qs the pair creation efficiency of the SPDC crystal 
and At equivalent to the detection gate width of the 
heralded photons detection. We can now express F by 

P{l)q H VH 
F - At— 

which allows us to see which parameters should possibly 
be improved. We note that F does not depend on the 
efficiency of the SPDC crystal. Thus, using more efficient 
pair generation, e.g. by using waveguides in periodically 
poled crystals |l5|. would not impact on the quality of 
this single photons source. However, the source can be 
improved by: i) Increasing P(l) by using better optics, 
coatings or alignments |l4| . or using a different photon 
pair creation technique as the Kerr effect inside optical 
fiber |2l| : ii) Improving the coupling and detection effi- 
ciency of the heralding detector; iii) Reducing the timing 
jitter of the heralding detector which would allow to use 
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much shorter detection gate width At for the heralded 
photons. 

We found interesting to compare in a standard vi- 
sual way the antibunching properties of the non-classical 
light coming from our source to other published results. 
Therefore, we devise a way to use the standard technique 
used in the case of synchronous single photon source, even 
though our source is asynchronous. For synchronous sin- 
gle photon sources, the standard way is to use detection 
signals coming from the Hanbury-Brown & Twiss setup 
as start and stop for a time to digital converter (TDC) 
device. Indeed, the time between successive possible cre- 
ations and detections of photons is constant (fixed by a 
given clock signal), and such setups enable one to directly 
obtain histograms where bin n represent the probability 
of having two detections in different detectors, that hap- 
pened at an interval of n triggers. For true single photon 
source, the central bin (n = 0) is of course zero while 
the other bins behave as for a poissanian source. The 
deepness of this central dip allow one to have a intuitive 
information about the quality of the single photon source 
considered. 

For the present source, single photons are created at 
random times because of the continuous pumping and 
spontaneous non-linear process used. The elapsed time 
between a detection in one detector and the other is 
therefore completely random, and for this reason it is not 
possible to use a TDC. Instead, the number of triggers 
between detection in different detectors is computed di- 
rectly from the raw binary data by a home made software. 
Our reduction algorithm enables us to present standard 
histograms. 

More precisely, the algorithm is the following: the pro- 
gram sequentially reads the raw data and registers in bins 
how many times a start and a stop were separated by n 
triggers (this number is called M(n) in the following). 
As a convention, n is positive if the start was given by 
detector A and the stop by detector B, and negative in 
the opposite case. The bin n = corresponds to the 
case for which two photons were detected simultaneously 
(i.e. for the same trigger). As for standard TDC, there 
could be some invalid start (another start for the same 
detector before a stop on the other) that will be ignored. 
Moreover the search for a stop is done only within a lim- 
ited range: if no stop is found before a certain number 
of triggers (we used 100) the search is cancelled and a 
new start has to be found in the subsequent data. All 
the processes leading to the histogram construction are 
summarized on figure 0] Note that this technique could 
be possibly applied to synchronous sources too. 

The theoretical form for M(n), when n ^ is given by 



M(n)=Cp A (l-p B )^p B 



(2) 



for n > and with pa and pb swapped for n < 0, where 
C is a normalisation constant which depends on the to- 
tal number of triggers, on the maximal n allowed, and 
on pa and ps- Detections that do not occur simultane- 
ously are indeed coming from different photons pairs and 
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FIG. 4: Samples of raw data. Here, as an example, data 
reduction leads to a start and a stop separated by 4 triggers. 



thus lead to independent detection events at detectors 
A and B. The probability p(n) of having n triggers be- 
tween a start and a stop is thus given by p(n) = M ^ 1 ' , 
C is obtained by fitting the experimental results with 
the formula above for n ^ 0. Figure [S] shows a typical 
result for the whole histogram created by the algorithm 
described above, using 100 as range limit value. The cen- 
tral bin (corresponding to the measured value for M(0)) 
is nearly indicating a very low p(0) = ^ pab ■ The 
data are in a very good agreement with formula @ as 
shown by the fitted curve. The decreasing slope on both 
side of the center is related to the detection probabilities 
Pa and p B - Theses values are quite high in our experi- 
ment compared to others. Indeed, recall that pa and ps 
describe not only the detection efficiencies but also the 
collection efficiencies which are rather high in our case. 
Consequently, the slope in figure[S]is quite steep and may 
look unfamiliar. 
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FIG. 5: Typical resulting raw histogram displayed as points 
over the full range allowed by the algorithm. The plain curve 
is obtained by fitting the experimental values with the formula 
© and varying C . 



Another mean of presenting these data is to restrict the 
display to a limited bin range around the center, and to 
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normalize the values in such a way that 5» (2) (0) is directly 
readable on the graph. This is done by dividing all exper- 
imental values by C ' paPb- The height of the central bin 
corresponds then directly to g^ 2 '(0) without dark counts 
substruction (see appendix for details on dark count sub- 
straction) . Figure El shows this representation for the 
same data as on figure El 
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FIG. 6: Typical resulting normalized histogram, restricted 
to small values around the center. It is possible to infer the 
(0) values from the central bin. This histogram correspond 
to acquisition made for 1.6 mW of pump. Note that the bins 
are made using raw values (dark counts are not substracted) 
and that in this case g (2) (0) = 3 x 1CT 3 . 



Note that it is in principle possible to turn our source 
into a synchronous one. A first step would be to use a 
pulsed laser to pump the non- linear crystal. A second 
step, following |23, would consists in a switched storage 
delay line (see figure 0). The latter would allow one to 
increase the probability of getting precisely one photon 
provided that the switch has a low enough insertion loss. 
The heralding signal is activated each time the storage 
loop is released (e.g. once every 10 pump pulse cycles), 
provided the Si-APD did register one photon. 
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FIG. 7: Schematic for using a SPDC source for high efficiency 
synchronous single photons production 



We demonstrated that a simple and well known SPDC 
setup 0, can be used to implement a very efficient 
high quality single photon source. The properties of 
our source for different pump powers have been inves- 
tigated, showing that good tradeoffs between very low 



<7 ( - 2 - ) (0) and high enough repetition rate can be found, 
which are competitive with much more complex sources. 
Moreover, SPDC sources can be designed to produces 
photons at various wavelength to match particular needs. 
Our source produces telecom wavelength single photons 
that are already in a standard single mode telecom fiber, 
is simple to operate and works at room temperature with- 
out stabilization. It is therefore a very good candidate for 
practical real-world telecom applications for which single 
photons are required. 



APPENDIX 

Let us elaborate on the way we obtain the relevant 
source parameters after substraction of the dark counts 
for our setup (see figure 
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FIG. 8: Schematic of the quantities used for calculation of the 
source parameters. 



The probability p count for a single detector to register 
a count is given by 

p cou„t = ^ _ p dcj p phot + p dc (! _ p phot) 

= p phot (1 - 2p dc ) + p dc 

where p dc is the probability for the detector to regis- 
ter a dark count in the absence of photons, and with 



P 



phot 



J2j>o Pj (1 — 0- ~ 7 ?) J ) t ne probability of detect- 
ing photons. Here Pj Pj = 1) represents the photon 
numbers distribution directly at the input of the detec- 
tor which features a quantum efficiency r\. At the input 
of the test bench, the beam splitter randomly separate i 
incoming photons in the two arms. Pj is thus the bino- 
mial probability that, having i photons before the perfect 
beam splitter, j photons reach a given detector (j < i) 

i.e. Pj = Q^ra^r = 2 ~ l Q- Usin S this P i into the 
above formula for pP hot , and using Q) 2 ^ = (1 + X Y> 
p phot reduces to p phot = 1 - (1 - | )\ Therefore, when- 
ever a heralding signal announces the presence of a multi- 
photon state (with probability of having i photons given 
by P(i), P(i) = 1) at then input of the beam splitter, 
Pa (the expression for pb is similar) is expressed by 



PA 




2p 



A ) 




(A.1) 



Here 77 is substituted by (3a, the overall efficiency detec- 
tion including fiber losses and quantum efficiency. 
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In the same manner, the joint probability pab can be 
calculated to be 



i=o 



p AB = £ P(i) £ ( * j -{ [(1 - (1 - PaY) (1 - 2p A C ) 
i=0 j=0 

x [(l-(l-f3 B f-rt)(l-2pi c )+pi c ] } 

=E^){( 1 -( 1 

/3a + /?s 



■4 C 



Pa 
2 




(1 - 2p A c ) (1 - 2p| c ) E2 



(1-2PA C )PB 

{l-2p% c )p A c 



CE2fc) 



The lines <|A.2fa,(l correspond to the cases for which at least 
one photon is detected simultaneously in each detectors, 
and are consequently non-zero only for i > 1. The lines 
(IA.2b|l are the terms corresponding to the cases for which 
at least one photon is detected by one detector while the 
other one registers a dark count. These expression are 
therefore non-zero only for i > 0. The last line (|A.2b|) 
corresponds to the case for which both detectors register 
dark counts, and is of course independent of i. 

Assuming P(0) + P(l) + P(2) = 1 in (jA~T|) and lfO|) 
(which is similar to neglecting terms for i > 2), and solv- 
ing for P(l) andP(2) (P(0) = l-P(l) -P(2)), it is pos- 
sible to get values for these quantities from a function of 
Pai Pab> Pa/b an d P A a/b (™ distinguish measured values 
for a parameters). P(l) and P(2) were also calculated 
from pg and p™ B and the final published values are the 



averages. These values take directly into account the ef- 
ficiencies and the dark counts and are thus as close as 
possible to the AHSPS parameters. 

The net values for pa, Pb and pab, denoted with nct , 
are obtained from 



Pa 



Pa 



Pa 



Tl I1Ct 

Pab 



Pab 



Pa-Pe 



by evaluating the right members of the equalities using 
the value for P(l) and P(2) calculated from p A , 
Pa/b and P d A/B ( or 'Pb instead of for p A et )- These 
values are used to calculate the net g^> (0) of table |U 
using ^ 

A simplified case with P% <C Pi and with p A c = pj^r = 
(no noise) and Pa = Pb = (3 (these conditions correspond 
to a good approximation to our experimental results), 
leads to: 



and thus 



iff =1^=^(1)1 , PAB=pmY 



P(1) = 2 P -A = PA±Eb > P{2 ) = 2 P -p- 



In this case we have 



( 2) ^_ Pab _ Pab = 2 P ( 2 ) 



P'XP'B 



(P 



m\2 



p(iy 



and thus we obtain the well known expression for (0) 
which enable one to compare the present source with a 
Poissonian source (as c/ (2) (0) = 1 for such light sources ) . 
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